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ABSTRACT 

The three-dimensional incompressible, steady and laminar flow field 
around a prolate spheroid at incidence Is considered. The parabolized Navier- 
Stokes equations are solved numerically. The method can handle vortex type as 
well as bubble type flow separation because the pressure is one of the 
dependent variables. In the present paper the distribution of the skin 
friction Is reported for two test cases. The first test case is a prolate 
spheroid of aspect ratio of 4:1 at 6° incidence and Reynolds number of 10^ 
(based on half the major axis). The second case is a spheroid with a 6:1 
aspect ratio at 10° incidence and Reynolds number of 0.8» 10^. The properties 
of the flow field near the body are discussed on the basis of the pattern of 
the skin friction lines, and the shape of the separation lines. Favorable 
agreement with experimental results is obtained. 
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Introduction 


The analysis of incompressible three-dimensional flows around slender 
bodies at incidence is very important in modern aerodynamics and 
hydrodynamics. The flow field is characterized by crossflow reversal, 
localized thickening of the boundary layer, formation of longitudinal vortices 
regions of backward flow and strong viscous-inviscid interaction*. The 
complicated phenomenon of flow separation is of special importance to 
aerodynamic designers as it has a major influence on the aerodynamic 
coefficients. In two-dimensional flows the point of separation from a surface 
coincides with the point at which the wall shear stress vanishes. In three- 
dimensional flows, such a criterion does not exist. Indeed, the shear stress 
does not vanish on the separation line except, possibly, at isolated points. 

At present, the understanding of the three-dimensional flow separation relies 
mainly on observations drawn from the analysis of the flow pattern near the 
body. Particularly fruitful results were obtained by visualization of the 
pattern of the skin friction lines on the surface of wind tunnel models with 

O 

oil streak techniques . The theory of the skin friction lines is based on 
Lighthill's work"* which shows that the skin friction lines have to obey 
certain topological rules. It is widely accepted now that a necessary 
condition for flow separation in three-dimensional flow fields is the 
convergence of many skin friction lines into a single line*. Whether it is 
also a sufficient condition is a matter of current debate. In a similar 
manner, the divergence of the skin friction lines from a single line 
characterizes the attachment or reattachement of the flow. 

We have already mentioned the difference between the nature of separation 
in two-dimensional and three-dimensional flows. This has given rise to 
various definitions of separation. Thus Masked used experimental evidence 
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to classify two types of flow separation in three-dimensional flow fields: 
bubble type separation and free vortex type separation. In the bubble type 
separation the main direction of the flow is reversed in the separated zone. 

In the vortex type separation no velocity reversal occurs in the streamwise 
direction. The flow leaves an oblique separation line, usually as a free 
shear layer, and rolls up along the downstream direction. Thus only 
circumferentially reversed flow exists. Wang^ defines two kinds of flow 
separations as well, namely "closed" separation and "open" separation. These 
terms are based on the shape of the separation line. The closed and open 
separations are very similar to Maskel's bubble and vortex type separations, 
respectively. Other investigators prefer to describe the bubble type 
separation as "two-dimensional" separation, indicating the reversal of the 
mainstream direction, which is typical of separated regions of two-dimensional 
flow fields. The vortex type separation is called "three-dimensional" 
separation to emphasize the fact that this kind of flow separation may occur 
only in three-dimensional flow fields. 

As the phenomenon considered is so complicated much of the research was 
restricted to simple geometries. One of the simplest possible geometries is 
the prolate spheroid, which has been studied experimental ly® - ^ and 
numerically® - ^. In all the numerical calculations the boundary layer 

O 

approximation was used. For instance, Wang solved the degenerate case of 
the boundary layer equations at the windward and at the leeward symmetry 

Q 

planes. Later, Wang computed the flow field around a spheroid of major to 
minor axis ratio of 4:1 at 6° incidence. A solution was obtained over a 
significant part of the body. At the rear part of the spheroid no solution 
could be found because of numerical instabilities. Wang concluded that the 
flow is separated there. Cebeci et al.^ 0 obtained very similar results for 
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the same case and calculated the flow at other angles of attack as well. In 

the above computations, the boundary layer equations were solved by marching 

along both the axial and circumferential directions, but Hirsch and Cebeci^ 

12 

and Patel and Choi solved the momentum equations by an ADI method which 
was expected to be more stable in circumferentially reversed flow regions. In 
addition, Hirsch and Cebeci^ also solved what they call the 
"Parabolic-Elliptic Boundary Layer" equations in which the pressure field is 
known, but the circumferential diffusion terms are retained. No significant 
change in the results was reported. 

The boundary layer equations for a 6:1 spheroid at incidence were solved 
by Schoenauer^ and Patel and Baek^. The code written by Schoenauer^ 
is for general bodies and uses adaptive high order methods with control on the 
truncation error bounds. Unfortunately, no detailed results which can be 
easily compared with other results are given. Stock^ solved the integral 
boundary layers equations for a variety of cases. Not surprisingly, 
reasonable agreement with the experiment was obtained only for those regions 
of the flow field in which the flow direction does not deviate significantly 
from that of the main flow. 

In all the above mentioned boundary layer calculations flow fields without 
flow reversal were readily calculated. These methods may (with special 
treatment) handle also weak circumferential flow reversal, provided that 
separation does not occur. They can not handle separated flow of both kinds. 
In general the boundary layer approximation depends on specification of 
pressure, say by a potential flow solution.. In such calculations the 
separation line cannot be reached, and the formation of numerical difficulties 
is usually interpreted as the onset of separation. Against this background we 
have to note the success of the thin layer approximation^ in the 
calculation of compressible separated flow on slender bodies at incidence. 

This is attributed to the fact that the pressure is not prescribed. 
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A solution of the full Navier-Stokes equations is very expensive and has 
been attempted only by few investigators. Haase^ solved the full equations 
using vorticity /velocity formulation, and a cartesian coordinate system. As 
the coordinates were not body-fitted special techniques were used near the 
wall. His results appear to suffer from too low a resolution, and we believe 
that his technique requires the addition of a potential velocity correction to 
ensure mass conservation. 

In this paper, the laminar, steady and incompressible flow field around a 
prolate spheroid at incidence is investigated using numerical solutions of the 
parabolized Navier-Stokes equations in which the pressure is not 
predetermined. Consequently, the equations are not singular at separation 
lines and separated flow regions may be included in the solution domain. In 
the present work we consider the skin friction distribution only. The pattern 
of the skin friction lines is computed and particular attention is given to 
the separation lines. 

Mathematical Model and Numerical Solution 

The flow is approximated by the steady and incompressible Parabolized 
Navier-Stokes (PNS) equations which are obtained from the full Navier-Stokes 
equations by neglecting the streamwise diffusion terms in the momentum 
equations. This stands in contrast to the parabolic approximations to the 
Navier-Stokes equations (or the boundary layer equations) in which the 
pressure is given a priori. As the pressure is not predetermined in the 
parabolized Navier-Stokes equations the problem of singularity near separation 
(which is typical to the boundary layer equations) does not arise here. 

The governing equations are formulated in a general axisymmetric 
curvilinear orthogonal coordinate system, using primitive variables, and 
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contravariant velocity components. The three coordinates p, e, 5 (see Fig. 1) 
run in the normal, circumferential and the main stream direction, 
approximately. These coordinates are stretched to q, s, t coordinates 
respectively which are orthogonal as well. The equations are simplified 
neglecting the stream-wise diffusion terms. The equations are: 
the continuity equation: 
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where 


M = oU f^ + xV fl + 6w ft 

( 2 a) 

a = — 5 — a = a J 

h . 

( 2 b) 

6 = x ** BJ 

h s 

( 2 c) 

Y = 6 = yJ 

h t 

( 2 d) 

and J is the Jacobian: 


J = h„ h„ h* 
q s t 

( 2 e) 

The Lamme coefficients are given by: 


9 3 aX. 2 

h Q = fLi 

(3a) 

, 3 sX, 2 
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(3b) 
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where (X^, X£, Xg) is the cartesian coordinate system. 


The contravariant velocity components V^, V s , V^. were 

transformed to: 

U " h q V Q 

(4a) 

V = h s V $ 

(4b) 

« = h t V t 

(4c) 
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In the present study, a prolate spheroidal coordinate system was chosen 
with the following Lamme coefficients: 


do(q) 

dq 

sh 2 p(q) + sin 2 e(t) 

(5a) 

de(s) 

ds 

shp(q)sins(t) 

(5b) 

dt 

sh 2 p(q) + sin 2 j-(t) 

(5c) 


p(q), e-(s), (t) are one-dimensional stretching functions. 

The computed flow field did not cover the entire prolate spheroid. The 
upstream boundary was placed some distance downstream of the forward 
stagnation point at a region where the boundary layer approximation can still 
be justified and can provide the upstream boundary conditions for the PNS 
equations: 


u = u 


up’ 


v = v 


up’ 


w = w 


up 


(6a) 


where the subscript "up" stands for upstream conditions. 

v 

The downstream boundary was placed ahead of the rear stagnation point in 
order to minimize the usage of computational resources. It should be 
emphasized that the PNS equations are not signular at reversed flow regions, 
yet the bubble type separation at the rear part can be accurately computed 
only if the downstream boundary is moved far enough into the wake. Moreover, 
the calculation in the separated region was found to increase the number of 
interations. Thus the total demand of computer resources became higher than 
our computer (IBM 3081D) could cope with. Due to the parabolization it is 
necessary to specify only one condition at the downstream boundary. We chose 
to specify a pressure boundary condition as follows: 


IE = (IE) 

at 'at ; pot 


(6b) 


where "pot" stands for the potential value. 
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The outer boundary was placed in the potential flow region. Three 
boundary conditions must be specified there, for example: 

v " v pof w = w pot> P = Ppot < 6c > 

On the body the usual no-slip and no injection conditions were applied: 

u = v = w = 0 (6d) 

Because of the symmetry of the flow, the solution was obtained for a half of 
the field only. On the leeward and windward symmetry planes symmetry 
conditions were used. 

The finite difference equations written over a staggered grid were solved 
iteratively by a consistent and stable procedure without any further 
approximations. More details on the numerical method were reported by 
Rosenfeld and Israeli* 8 for the two-dimensional case and by Rosenfeld*^ 
for the three-dimensional case. 

A typical grid consisted of about 15000-20000 points: 25 in the normal and 
circumferential direction and 25 to 33 points along the spheroid. About one 
to two hours of CPU time were needed for a solution on an IBM 3081D computer, 
depending on the number of grid points and on the angle of attack. No 
standard convergence tests could be performed by systematically increasing the 
number of grid points because of the excessive computational cost. 
Nevertheless, some partial accuracy tests were conducted by varying the 
location of the outer and downstream boundaries while keeping the total number 
of mesh points constant. No significant differences were recorded for the 
cases reported in this paper. 
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Results 

The flow field was solved about two slender prolate spheroids of major to 
minor axes ratio of 4:1 and 6:1 and for several angles of attack. For the 
sake of brevity, the results of only two cases will be reported here, one for 
each thickness ratio. All the calculations are for laminar flow. The 
Reynolds numbers chosen are such that laminar flow may be expected to prevail 
over much of the flow field (except at separated zones). In the presentation 
the streamwise coordinate is the normalized axial distance along the spheroid 
z (see Fig. 1): z=-l and z=l correspond to the nose and to the rear end of the 
spheroid respectively. The angles e=0* and ©=180* correspond to the windward 
and leeward symmetry planes respectively. 

3.1 A 4:1 spheroid at 6* incidence 

Several numerical solutions of this case using the boundary layer 

Q 1 O 1C 

equations were reported in the literature ’ . Unfortunately, no experi- 

mental results are available for comparison with the calculations. Yet, 
comparison with previous boundary layer solutions is useful because at such a 
low incidence, the extent of the separated regions is very limited and 
therefore the boundary layer solutions may be expected to be quite similar 
to the parabolized Navier-Stokes results over a significant part of the 
spheroid. The Reynolds number based on half the length of the major axis of 
the spheroid and on the uniform velocity upstream of the body was set to 
10®. The upstream boundary was placed at z=-0.8 and the initial conditions 

O 

were approximated from the boundary layer solution of Wang for the same 
case. 

The distributions of the skin friction coefficients are compared to the 
boundary layer approximation results of Wang® and Patel and Choi^ in 
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Figs. 2 and 3. The skin friction is normalized by the Reynolds number: 

C F = /Re* (7) 

•jpu 

where Re is the Reynolds number and T w is the shear stress on the wall, p is 
the density and U is the free stream velocity. Figures 2a and 2b show the 
distributions of the skin friction coefficient on the windward and leeward 
symmetry planes, respectively. Examination of the results show that the two 
boundary layer solutions are not identical. The difference between the two 
boundary layer solutions in both the windward and the initial part of the 
leeward side is fairly uniform. Our results lie generally between the two 
boundary layer solutions, although we used Wang's solution as the initial 
condition. However, in the downstream part of the leeward side the two 
boundary layer solutions approach one another, whereas the present value of 
the skin friction is somewhat higher. This region is characterized by the 
thickening of the boundary layer and ultimately separation occurs. The 
interaction between the separated flow and the pressure field is expected to 
lead to quite substantial departures from the potential pressure field 
specified in the boundary layer calculations, and consequently the 
disagreement between the present and boundary-layer results might have been 
expected. 

In the present work, the separation point in the leeward plane is found to 
be at z=0.75 (see Fig. 2b). The separation points computed by the boundary 
layer approximation are quite close to one another: z=0.724, z=0.72 and z=0.73 
for Wang®, Cebeci et al.*® and Patel and Choi^ respectively. Once 
more, the difference may be attributed to the viscous-inviscid interaction 
which is neglected in the boundary layer approximation. 
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The circumferential distribution of the t and s components of the skin 
friction coefficient are shown in Fig. 3 at three axial locations, and 
compared with the data of Patel and Choi^. As the agreement of the present 
skin friction coefficient on the planes of symmetry with Patel and Choi is not 

very good (see Fig. 2) it is not surprising that the present axial component 

of the skin friction in the frontal parts of the spheroid differs from the 
boundary layer one by a fairly uniform value. Another source for the 
disagreement can be found in the coarse circumferential resolution used by 

Patel and Choi: intervals of 20* in contrast to 7.5* in the present 

calculation. Still the circumferential component of the skin friction is in 
good agreement. The situation is different in the rear part of the spheroid 
where circumferentially reversed flow prevails and axial separation appears to 
be imminent. 

The distribution of the calculated skin friction coefficient vectors on 
the unwrapped surface of the spheroid is shown in Fig. 4. Three distinct 
regions can be observed. In the first region, the vectors point towards the 
back and the leeward side of the spheroid. Further downstream and on the 
leeward side the direction of the vectors are shifted towards the windward 
side, suggesting the onset of a vortex. At the very end of the spheroid, the 
vectors point towards the nose, indicating the formation of a separation 
bubble. The solid line marks the location where the circumferential component 
of the shear stress reverses its direction. 

The skin friction lines are tangent to the shear stress vector on the 
surface and can be computed^ - from: 

h s ds _ C F,s 

h t dt C F,t 


( 8 ) 
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As the skin friction lines are the projections of the streak lines near the 
body*, their pattern may be instructive of the flow field in the vicinity of 
the body. Lighthill has shown that in the three-dimensional case the flow 
may separate from a solid body because of two reasons: Cp ->0 and/or the 
convergence of the skin friction lines into one particular line. An 
attachment (or reattachment) line of the flow is characterized by the 
divergence of the skin friction lines away from the attachment line. 

The skin friction lines were computed in the present work by solving (8) 
with a second order Runge-Kutta method. The components of the skin friction 
are interpolated using a bilinear spline approximation. The skin friction 
lines for the 4:1 spheroid at 6° incidence are shown in Figs. 5. Figure 5a 
shows the pattern on the unwrapped spheroid while the side, bottom and upper 
views of the skin friction lines on the spheroid are shown in Fig. 5b. Two 
types of skin friction lines convergence are evident: a very short swept line 
and a second line that appears to encircle the spheroid at about z=0.70. The 
second line of convergence is obviously a separation line that divides the 
flow into two regions: a region that is accessible to the flow from the 
forward stagnation point and a region that is not accessible and is in fact a 

Q 

bubble type separated region with reversed flow. In the terms of Wang this 
separation is called "closed" as it originates from a closed separation line. 
Figure 6 compares the zero circumferential shear stress line and the closed 

q IQ 

separation line with the results of Wang and Cebeci et al. . The 
qualitative shape of the closed separation line is similar to other 
computations although the "tounge" of the separation line in uur results is 
less pronounced and more close to the windward side. It should be noted that 
all separation phenomena are not likely to be well predicted by the boundary 
layer approximation. The agreement of the zero circumferential shear stress 
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line is much better, indicating that the boundary layer approximation is valid 
in regions not close to flow separation. Although there is circumferential 
flow reversal from about z=-0.1, the skin friction lines converge only at 
z=0.45 about 15* further in the leeward direction to the zero circumferential 
shear stress line (see Figs. 5a). However, the pattern of the skin friction 
lines is not sufficient to conclude whether a vortex type separation has 
occurred there. 

The topology of the skin friction lines found in our calculations is 
notably similar to that conceived by Han and Patel® from their flow field 
visualizations at 10° incidence and Reynolds number of 4.10 4 (see Fig. 10 in 
their paper). Cebeci et al.^ suggest three possibilities of the shape of 
the swept convergence line near the closed separation line. The present 
results are similar to their second option (Fig. 10 in their paper). 

3.2 A 6:1 spheroid at 10° incidence 

In this higher incidence case a bubble type separation exists as well as a 
free-vortex type separation. Some boundary layer solutions were reported for 
this case*® - *®. The experimental data available is reported by Kreplin et 
al.^ who measured the skin friction shear stress on a 6:1 spheroid at a 
Reynolds number of 8*10®. The same Reynolds number was used in the present 
calculations. 

The absolute magnitude of the normalized shear stress in several z 
stations is compared in Fig. 7 to the experimental values^ and to the 
computed boundary layer results of Patel and Baek*^. The agreement at small 
distances from the upstream boundary is not good, possibly due to poor 
specification of the initial conditions at the upstream boundary in the 
present calculations which were approximated from Stock's integral boundary 
layer solution*®. Further downstream the agreement of our calculation with 
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the experimental results improves while the results of Han and Baek get 
worse. Kreplin et al. related the high values of the skin friction near the 
leeward side to the transition of the separated flow to turbulent flow. In 
our laminar calculation this phenomenon cannot be reproduced. 

Figure 8 shows the shear stress vector plot on the unwrapped spheroid. A 
minimum along the circumferential direction in the shear stress is clearly 
seen. It is interesting to note that due to three dimensional effects, this 
minimum does not occur along the leeward side where the potential adverse 
pressure gradient is maximal. Figure 9a shows the shear stress vector plot 
interpolated from our computations with the experimental results of Kreplin et 
al.^. In the laminar region of the flow field the agreement is usually good 
bearing in mind the complexity of the flow field and the difficulties in 
measuring the stress in the laminar regions as reported by Kreplin et al. 
However, the agreement at the initial stations near the leeward side is not as 
good. This is attributed to the poor initial conditions used in the present 
calculation at the upstream boundary. Figure 9b shows the same comparison 
with the boundary layer results of Schoenauer^. Boundary layer solutions 
were not obtained for significant parts of the flow field, although these are 
the most interesting regions where separation may occur. In the other 
regions, the agreement with our results is better than with the experimental 
results. 

Figures 10 show the skin friction lines on the surface. In this case no 
bubble type separation is found in the region solved. On the other hand, a 
long and swept convergence line is present. The precise origin of the line 
cannot be defined with certainty. The skin friction lines merge into a line 
of convergence from both sides. It is interesting to note that now the zero 
circumferential shear stress line is only a small distance in the windward 
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direction from the convergence line, /'gain, it is impossible to predict the 
onset of an open type separation based on the properties of the flow field in 
the vicinity of the surface. 

Kreplin et al. obtained the skin friction lines from the measured 
distribution of the shear stresses (Fig. 5 in [7]). The main separation line 
appears to be placed further in the leeward direction than the line we have 
computed. The differences may be explained by the transition to turbulent 
flow that occurred in the experiment, as turbulent flow can withstand higher 
adverse circumferential pressure drop before separation. Kreplin et al.^ 
observed a second separation line very close to the leeward side. We have not 
found this line, although there is a tendency of creating a second convergence 
line in the rear part of the spheroid. 

Conclusions 

The incompressible, steady and laminar parabolized Navier-Stokes equations 
were used to simulate the flow field around a prolate spheroid. Solutions 
were obtained not only in regions of attached flow, but also in domains where 
the flow is separated in the lateral and/or the axial direction. In this 
paper we concentrate only on the skin friction on the surface of the 
spheroid. Even so, this information is useful as it gives clues on the nature 
of the flow, and allows comparisons with available experimental and numerical 
data (anyhow, much of this data is confined to the skin friction). The 
results obtained agree with previous experimental data and numerical solutions 
of the boundary layer approximation, in the non-separated regions.’ In the 
separated region the skin friction is too low compared to the experimental 
data, and this is attributed to the turbulent flow prevailing there. Yet, the 
capability of the method to yield prediction in the separated regions has been 
demonstrated. 
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The present calculations can be used to study the conditions and symptoms 
of separation. Further analysis of the data in the whole flow field by 
appropriate visualization of the numerical solution is needed to obtain final 
conclusions. -It should be noted that complete simulations of 
three-dimensional separated flow fields require turbulent flow modelling as 
well, since separated regions are usually turbulent. Also, the boundary 
conditions at the upstream boundary must be improved, perhaps by solving the 
flow field from the forward stagnation point. These studies are currently 
under way and will be reported in due course. 
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